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*~i ' Abstract. We introduce a metric notion of Ricci curvature for PL 

manifolds and study its convergence properties. We also prove a fitting 
version of the Bonnet-Myers Theorem, for surfaces as well as for a large 
class of higher dimensional manifolds. 



1. Introduction 



q 

r~| ' Recent years have seen a great "revival" of Ricci curvature, due mainly to 

C^ ' Perelman's celebrated work on the the Ricci flow and the Poincare conjecture 

[Per02j ■ |Per03| . but also to its extension to a far larger class of geometric 

objects, than merely smooth (3-)manifolds (see |Vi09j and the bibliography 

►-j^ therein). In consequence, Ricci curvature has become an object of interest 

CN ' and study in Graphics and Imaging. The approaches range from implemen- 

0\ ' 

ly-^ ' tations of the Combinatorial Ricci curvature of Chow and Luo |CL03| - see, 

e.g. [GY08] ■ through classical approximation methods of smooth differen- 
tial operators |AslO| . [FAF09] . to discrete, purely combinatorial methods 
Csl ; [SAWZOQj . 

We have addressed the problem of Ricci curvature of PL surfaces and 
higher dimensional PL (piecewise flat) manifolds, from a metric point of 
view, both as a tool in studying the Combinatorial Ricci flow on surfaces 
[Sallaj . and, in a more general context, in the approximation in secant of 
curvature measures of manifolds [Sallbj . Computational applications aside, 
these and related problems - see |Ru07] . [BernOSa] - make the study of a 
robust notion of Ricci curvature for PL spaces a subject of thriving interest 
in the Geometry and Topology of (mainly 3-dimensional) manifolds. 



Date: February 25, 2013. 

1991 Mathematics Subject Classification. 51K10, 53C21, 65D18. 

Key words and phrases. PL manifold, Ricci curvature, Bonnet-Myers Theorem, Wald- 
Berestovskii curvature. 

Research supported by European Research Council under the European Community's 
Seventh Framework Programme (FP7/2007-2013) / ERG grant agreement n° [203134]. 

1 



2 EMIL SAUCAN 

This article represents a continuation of (or, in a sense, an appendix to) 
both papers above: In the first part of this paper we address the main prob- 
lem, namely that of defining a computable, discrete metric Ricci curvature 
for PL (piecewise flat) manifolds. Here we use methods similar to those in 
[Sallbj . to address a problem - or, rather, a particular case - considered 
therein. We also investigate the convergence properties of this newly intro- 
duced curvature. In the second (and last) part we address a rather more 
theoretical question, namely if the newly introduced version of Ricci curva- 
ture satisfies - as indeed expected from a proper ( "correct" ) notion of Ricci 
curvature - a Bonnet-Myers type of theorem. (The methods in this part 
are, partially, those developed in |Salla| .) 

A note to the reader before we proceed to the main part of our paper: 
Our default source for geometric differential definitions and results is |Be03j , 
and, if no other source is specified the reader should consult, if needed, this 
encyclopedic source. Also, as background material for PL topology we refer 
the reader to [Hu69j . 

2. Definition and Convergence 

To begin with, we have to be able to properly define Ricci curvature for 
PL manifolds. This is indeed possible, not just for PL manifolds but also for 
polyhedral ones - and in a quite natural manner - combining ideas of Stone 
|St76aj . |St76b] and metric curvatures. For this one regards Ricci curvature 
as the mean of sectional curvatures: 

n 

(1) Ric(ei) = Ric(ei,ei) = '^K{ei,ei), 

i=2 

for any orthonormal basis {ei,--- ,e„}, and where K(ei,ej) denotes the 
sectional curvature of the 2-sections containing the directions ei. 

First, one has, of course, to be able to define (variational) Jacohi fields. 
This is where we rely upon Stones's work. However, we do not need the 
whole force of this technical apparatus, only to determine the relevant two 
sections and, of course, to decide what a direction at a vertex of a PL 
manifold is. 

In fact, in Stone's work, combinatorial Ricci curvature is defined both 
for the given simplicial complex T, and also for its dual complex T* . In 
the later case, cells - playing here the role of the planes in the classical 
setting of which sectional curvatures are to be averaged ~ are considered. 
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However, his approach for the given complex, where one computes the Ricci 
curvature Ric(o", ri — T2) of an n-simplex a in the direction of two adjacent 
(n — l)-faces, ti,T2, is not natural in a geometric context (even if useful 
in his purely combinatorial one), except for the 2-dimensional case, where 
it coincides with the notion of Ricci curvature in a direction (i.e., in this 
case, an edge - see also Remark 12.81 below). Passing to the dual complex 
will not restrict us, since (T*)* = T and, moreover - and more importantly 
- considering thick triangulations enables us to compute the more natural 
metric curvature for the dual complex and use the fact that the dual of a 
thick triangulation is thick, as we shall detail below. Working only with 
thick triangulations does not restrict us, however, at least in dimension < 4, 
since any triangulation admits a "thickening" - see |Sa05] l I 
First, let us recall the definition of thick triangulations: 

Definition 2.1. Let r C M" ; < /c < n be a /c-dimensional simplex. The 

thickness (or fatness) 93 of r is defined as being: 

dist(6,(9a) 

(2) V{r) = — , 

diamo" 

where h denotes the barycenter of a and da represents the standard notation 

for the boundary of a (i.e the union of the (n — l)-dimensional faces of a). 

A simplex r is ipQ-thick, for some (pQ > 0, if (/^(r) > ipo. A triangulation (of 

a submanifold of M") T = {o"j}jgi is ip^-thick if all its simplices are ipQ-thidi. 

A triangulation T = {crjjjgi is thick if there exists v'o ^ such that all its 

simplices are ip^-thick. 

Remark 2.2. This is Munkres' definition |Mu66| . For a discussion of other, 
equivalent definitions, their mutual interplay and relationship with cer- 
tain aspects of Differential Geometry (mainly curvature approximation) see 
|Sallb| . Note that this definition holds for more general simplices, not nec- 
essarily Euclidean ones. 

To be able to define and estimate the Ricci curvature of T and T* and the 
connection between them, we have to make appeal in an essential manner 



This holds, as already mentioned, for any PL manifold of dimension < 4, and in all 
dimensions for smoothable PL manifolds, as well for any manifold of class > C^ . Since the 
proof of the main result of Section 3, regarding manifolds of dimension higher than 3, holds 
only for manifolds admitting smoothings, restricting ourselves only to such manifolds does 
not represent anyhow a great hindrance. 
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to the fatness of the given complex. We begin by noting - by keeping in 
mind formula ([2]) - that, since the length of the edge /*•, dual to the edge kj 
common to the faces /j, fj equals ri + Vj, the first barycentric subdivisioro of 
a thick triangulation is thicko Since the definition of thickness also makes 
sense for for general cells (see [SAZ12] . Definition 3.1), we can summarize 
the discussion above as 

Lemma 2.3. The dual complex of a thick (simplicial) complex is thick. 

Moreover, we have the following (common) Gromov-Hausdorff conver- 
gence property: 

Lemma 2.4. Let T, T* be as above. Then 

(3) lim (r)= lim (T*) , 

5{r)^o <5(r)^o 

where 5{T),5{T*) denote the mesh ofT,T*, respectively. 

Remark 2.5. It is important to stress here the crucial role of the thickness 
of the triangulation, as far as geometry is concerned: Thickness ensures, 
by its definition, the fact that no degeneracy of the simplices occurs, hence 
no collapse and degeneracy of the metric can take place. Moreover, in its 
absence no uniform estimates for the edge lengths can be made, hence con- 
vergence of (dual) meshes and, as we shall see shortly, of their metric Ricci 
curvatures, can no be guaranteed. 

Returning to the definition of Ricci curvature for simplicial complexes: 
Given a vertex vo, in the dual of a n dimensional simplicial complex, a 
direction at vq is just an oriented edge ei = fo^i. Since, there exist precisely 
n 2-cells, ci, . . . , Cn , having ei as an edge and, moreover, these cells form part 
of n relevant variational (Jacobi) fields, the Ricci curvature at the vertex v, 
in the direction ei is simply 

n 

(4) Ric(^) = ^K(c,). 

4=1 



needed in the construction of the dual complex - see e.g. [Hu69] 

For planar triangulations, and also for higher dimensional complexes (embedded in 
some R'^), one can realize the dual complex (also in R^) by constructing the dual edges 
l*j orthogonal to the middle of the respective hj-s. To show the thickness of the dual 
simplices, one has also to make appeal to the characterization of thickness in terms of 
dihedral angles (Conditions (1.15) of |CMS84| '). 
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Remark 2.6. Observe that the index "i" in the definition (JH) above runs from 
1, and not from 2, as expected judging from the classical (smooth) setting. 
This is due to the fact that we defined Ricci curvature by passing to the 
dual complex, with its simple but demanding (so to say) combinatorics. For 
the implications of this fact, see Theorem 1 2 . 1 1 and Remark 12.121 below. 

Remark 2.7. Note that we followed |St76aj only in determining the varia- 
tional fields, but not in his definition of Ricci curvature. Indeed, he considers 
a direction at a vertex vq to be the union of two edges ei, 62 in the dual com- 
plex, where ei = (woi ^1)) ^2 = ('Wii '^2) and the direction is determined by the 
lexicographical order. Then (according to |St76aj ) the relevant variational 
field are given by the 2n distinct 2-cells Ci, . . . , C2n, containing the edges ei 
and 62, 2n — 1 of them containing one, but not both of them. Hence, the 
Ricci curvature at v in the direction 6162 is to be taken as the total defect 
of these 2n — 1 cells. This approach is necessary in the combinatorial case. 
However, it is more difficult than our approach and it would produce un- 
necessary complications in determining the relevant analogues of the (n — 1) 
2-sections of the classical, smooth case. Moreover, it is quite possible that, 
in any practical implementation, the advantages obtained by considering 
larger variational fields would be countermined by "noise" added by consid- 
ering such order 2 (or larger) neigbourhoods of the given vertex. However, 
computing Ricci curvature according to this scheme is still possible, using 
our metric approach (but see also the following Remark 12. 8p . 



Remark 2.8. It is still possible to compute Ricci curvature according, more- 
or-less, to Stone's ideas, at least for the 2-dimensional case. Indeed, accord- 
ing to |St76b] 

2n-l 

(5) Ric(cr, n - r2) = 8n - ^ {N{l3j) \ pj < n or pj < T2; dim/3j = n - 2} 

where N{j3j) denotes the number of n-simplices a, such that f3j < a. 

This definition of Ricci curvature is a combinatorial defect onqj. This is 
evident from its expression, but made more transparent by the 2-dimensional 
case: Indeed, in this case, the simplices Pj are 0-dimensional, i.e. vertices, 
and N{Pj) is just the number of 2-simplices having Pj as a common vertex. 



presumably inspired by the classical definition of Gauss curvature as the angular defect 
at a vertex — see, e.g. [HC-V]. 
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hence Ric((T, ti — T2) represents nothing but the total combinatorial defect 
at these 2n — 1 vertices. 

In consequence, using the approach of the original proof of Hilbert and 
Cohn-Vossen |HC-V52] . (and following methods well established in Graph- 
ics, etc.), we can consider, instead of the combinatorial defect, the angular 
defect of the cell Cj dual to the vertex f3j. This, of course, applies both for 
our way - as well as Stone's - of determining a direction. 

However, this approach to the definition of PL Ricci curvature is far 
less intuitive (and apparently has lesser geometric content, so to speak) in 
dimension > 3. This is the reason why, for our present study, we have made 
use of the dual complex. 

Remark 2.9. Note that, up to this point, we have not yet defined the sec- 
tional curvature K{c) of a cell c (see however the discussion below). Nev- 
ertheless, regardless of the specific definition employed, we obtain, quite 
trivially, the following generalization of the classical curvature bounds of 
Riemannian geometry (compare also with [Bern03b] . Theorem 1): 

Theorem 2.10 (Comparison theorem). Let M = Mpj^ be an n- dimensional 
PL manifold, such that Kw{M) > Kq > 0, i.e. K{c) > Kq, for any 2-cell 
of the dual manifold (cell complex) M* . Then 



(6) Kw = Kq^ Riciy = nKQ . 

Moreover 



(7) Kw = Kq ^ scaV ^ n{n + l)Ko , 

where scalvv' denotes the scalar metric curvature of M, defined as scalw{v) = 
Y^ K\y{c), the sum being taken over all the cells of M* incident to the vertex 
vofM*. 

Remark 2.11. Note that inequality d?]) can be formulated in the seemingly 
weaker form: 



(8) Ricw F nKo =^ scalw = n{n + l)KQ, 
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Remark 2.12. Note that in all the inequalities above, the dimension n ap- 
pears instead of n — 1 as in the smooth, Riemannian case (hence, for instance 
one has in ([7]), n(n + 1)Kq, instead of n(n — l)-fCcl_| as in the classical case). 
This is due to our definition ^ of Ricci (and scalar) curvature, via the dual 
complex of the given triangulation, hence imposing standard and sirnple 
combinatorics, at the price of allowing for only for such weaker bounds|j 

To determine - using solely metric considerations - the sectional curva- 
tures K{ci) of the cells Cj, we shall employ the so called {modified) Wald 
curvature Kw (K^). At this point, we have to remind the reader a number 
of definitions and results that, unfortunately, are perhaps (at least partly) 
forgotten. We begin with following basic 

Definition 2.13. Let {M,d) be a metric space, and let Q = {pi,...,P4} C 
M, together with the mutual distances: dij = dji = d{pi,pj); 1 < i,j < 4. 
The set Q together with the set of distances {djj}i<jj<4 is called a metric 
quadruple. 

Remark 2.14. One can define metric quadruples in a somewhat more ab- 
stract manner, that is without the aid of the ambient space: In this ap- 
proach, a metric quadruple is defined as a 4 point metric space; i.e. Q = 
({pi, ...,pi],{dij}), where the distances d^ verify the axioms for a metric. 

We next introduce some necessary notation: Let 5^ denote the complete, 
simply connected surface of constant Gauss curvature k, i.e. 5^ = M^, if 
K = 0; £•« = §2 if /^ > 0; and 5« = M^ ^ if ^ < q. Here S^ = SV 
denotes the sphere of radius R = 1/\/k; and S'k = "^r^- stands for the 
hyperbolic plane of curvature \/—k, as represented by the Poincare model 
of the plane disk of radius R = Xj \J —k . Using this notation we can next 
bring 

Definition 2.15. The embedding curvature k{Q) of the metric quadruple 
Q is defined to be the curvature k of the gauge surface 5^ into which Q can 
be isometrically embedded. 

We are now able to bring the definition of Wald curvature |Wa35j (or 
rather of its modification due to Berestovskii |Ber86] ) : 



but even ii n — 3!... 
^without affecting the analogue of the Bonnet-Myers Theorem - see Section 3 below. 
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Definition 2.16. Let {X,d) be a metric space. An open set U C X is 
called a region of curvature > k iff any metric quadruple can be isometrically 
embedded in Sm, for some m > ki\ 

Remark 2.11 . Evidently, one can consider the Wald-Berestovskii curvature 
at an accumulation point of a metric space, hence on a smooth surface, by 
considering limits of the curvatures of (nondegenerate) regions of diameter 
converging to 0. 

Before we proceed further, let us make a certain modification of the no- 
tation, in order to make it more uniform and more familiar to the reader 
working in classical Differential Geometry as well as in Graphics: Henceforth 
we shall denote by Kw the Wald curvature of a surface {PL or smooth) , by 
analogy to its classical (Gauss) curvature K. (Of course, K\y{p) will denote 
the Wald curvature of a point on the surface.) 

At this point the question that naturally rises is whether it is possible to 
actually compute Wald curvature and, if possible, in what manner? First 
of all, the first, basic step is to note that the role of the abstract open sets 
U in Definition 12.161 above is naturally played by the cells Cj. We can state 
this as a formal definition, for the record: 

Definition 2.18. Let c be a cell with vertex set V^ = {vi, . . . ,Vp}. The 
embedding curvature K{() of c is defined as: 

(9) K{c)= min K{vi,Vj,Vk,vi) . 

l<i<j<k<l<p 

It is certainly worthwhile to note that it is possible to actually compute 
the Wald curvature of each of these cells, using the following formula for the 
embedding curvature n{Q) of a metric quadruple Q: 



(10) 
k{Q) 








ifr(Q) = 0; 














K, 


K<0 


if det{cosh -^ 


-K. 


dij) = 


= 0; 








K, 


fi;>0 


if det{cos y/H 


di. 


■) and 


\/k- 


dij 


< 


IT 



and all the principal minors of order 3 are > 0; 



While is not needed in the remainder of the paper, we mention for the sake of com- 
pleteness, that a metric space {X, d) is said to have Wald-Berestovskii curvature > «; iff 
for any x G X is contained in a region U of curvature > k. 
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where dij = d{pi,pj),l < i,j < 4, and T{Q) = T{pi, . . . ,p4) denotes, the 
Cayley-Menger determinant: 



(11) Tipo,...,P3) 



"10 ^ "l3 



"30 "'31 ■ ■ ■ 



Remark 2.19. (1) For some first numerical results regarding the appli- 
cation of these formulas in a practical context, see |Sa04j . |SA09j . 
However, it should be noted that, apart from the Euclidean case, 
the equations involved are transcendental, and can not be solved, in 
general, using elementary methods. 
(2) We have also employed Wald curvature as a malleable tool in con- 
junction with Ricci curvature in a somewhat more theoretical context 
in [Sallaj . On a more abstract note, we should remark that, given 
its (metric) intrinsic nature, K]y "behaves well", so to speak, un- 
der Gromov-Hausdorff convergence (see |BBIOO| , |Gr99] and |Sa04| , 
|SA09j for some applications in Graphics, Imaging, etc.). Moreover 
since it (or, rather a somewhat modified version of it Kw identifies 
with Rinow curvature (see |B153| . [BM70J ). it allows us to view the 
whole problem of defining and computing Ricci for PL (polyhedral) 
manifolds, (and in particular its applications in Graphics, Regge 
calculus, etc.) in the larger context of Alexandrov spaces (see, e.g. 
[BBIOOj . [Gr99]). 

Remark 2.20. Obviously, one can use the same method as above to com- 
pute the Ricci curvature (of T*), according to Stone's original approach for 
determining a directions in cell complexes. 

To return to the main problem of this section: From the definitions and 
results above we obtain - first discretely, at finite scale bounded away from 
zero - then passing to the limit) the following result connecting between the 
Ricci curvatures of a simplicial (polyhedral) complex and its dual: 

Theorem 2.21. Let T,T* be as above. Then 

(12) lim Ric(cj) = lim C-Ric*(cr*), 

mesh(7")— ^0 inesh(T*)— S-O 
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where a £ T and where a* £ T* is (as suggested by the notation) the dual 
of a. 

Remark 2.22. This result is, admittedly, somewhat vague. However, to 
our defense, we can only underline the fact that the precise constant C is 
hard to determine. The thickness condition, that ensures a metric "quasi- 
regularity" of the triangulation, supplies us only with weak estimates. To 
obtain stronger ones, one should be able to control the regularity of the com- 
binatoric structure, as well. (This is evident, but it will become even clearer 
in the sequel.) It should be noted in this context that, at least in Graph- 
ics, mesh improvement techniques allow us to consider such "combinatorial 
almost regular" triangulations. 

From Lemma [2.4l the fact that Ric(u) is defined in a purely metric, intrin- 
sic manner and from the fact that intrinsic properties are preserved under 
Gromov-Hausdorff limits (see |Gr99| ) and from Theorem 12.211 above, we 
easily obtain: 

Theorem 2.23. Let M" be a (smooth) Riemannian manifold and let T be 
a thick triangulation of M"^ . Then 

(13) RicT" — ;■ Ci • RicAf", as mesh(T) — ?■ , 

where the convergence is the weak convergence (of measures). 

For related results, see |CMS84j . |Sallb| for the Lipschitz-Killing curva- 
tures, [BK82J . [Sal la] for discrete (combinatorial, respective metric) Gauss- 
ian curvature, and |Bern06] . for the Einstein measures. 

Remark 2.24. While the desired constant Ci is, of course, Ci = 1, and some 
first experimental results hint that, at least for certain "nice" triangulations, 
this is indeed the case, we can't guarantee a better result - see the remark 
following the preceding theorem. 

Remark 2.25. While we have adopted the Wald curvature as the metric 
curvature for surfaceqj of our choice, for reasons detailed above, it would 
be interesting to explore the capabilities - both theoretical and practical - 
as far as PL Differential Geometry is concerned, of other metric curvatures 



and the Finsler-Haantjes one as a metric alternative for computation of principal 
curvature 
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(see [Sallbj and the bibliography therein) and in particular of the Menger 
curvature measure: 

(14) /x(r) = ^lp{r) = Y, 4AT)idi^^Tf , 

for some p > 1, where km denotes the Menger curvature (of the simplex T). 

3. The Bonnet-Myers Theorem 

Having introduced a metric Ricci curvature for PL manifolds, one natu- 
rally wishes to verify that this represents, indeed, a proper notion of Ricci 
curvature, and not just an approximation of the classical notion. Accord- 
ing to the synthetic approach to Differential Geometry (see, e.g. |Gr99j . 
|Vi09j ) . a proper notion of Ricci curvature should satisfy adapted versions 
of the main, essential theorems that hold for the classical notions. Amongst 
such theorems the first and foremost is Myers' Theorem (see, e.g., |Be03j ). 
And, indeed, fitting versions for combinatorial cell complexes and weighted 
cell complexes were proven, respectively, by Stone |St76aj . |St76bj . and For- 
man |Fo03j . Moreover, the Bonnet part of the Bonnet-Myers theorem, that 
is the one appertaining to the sectional curvature, was also proven for PL 
manifolds, again by Stone - see |St76c] . [St73j . 

For the special - yet of main importance in applications (see |CL03j . 
|GY08] . [Sal la] ) - case of 2-dimensional manifolds, such a result is easy to 
prove, given the fact that Ricci and sectional curvature essentially coincide. 
More precisely, we can formulate the following theorem: 

Theorem 3.1 (Bonnet-Myers for PL 2-manifolds - Combinatorial). Let 
Mp^ be a complete, closed 2-dimensional PL manifold of without boundary, 
such that 

(i) There exists do > 0, such that mesh(Mp^) < ddj 

(ii) Kcomb{Mli) >Ko>0. 

Then Mp^ is compact and, moreover 
(15) 

diam(M|,) < I ^-^ ^^, ^o>(2-V2)vr; 

^^' - 1 4vr3(io/[(2vr - do)(4^A:o - k^y/^ else ; 



Here niesh(M|.j^) denotes the mesh of the 1-skeleton of Mpi^, i.e. the supremum of 
the edge lengths. 
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r2 



where Kcomb denotes the combinatorial Gauss curvature of Mpj^, 

rap 

(16) Kcombip) = 27r - ^ ai{p) 

i=l 

where ai, . . . , amp are the (interior) face angles adjacent to the vertex Vi. 

Remark 3.2. Condition (IV that ensures that the set of vertices (of the PL 
manifold) is "fairly dense'Lj is nothing but the necessary and quite common 
density condition for good approximation both of distances and of curvature 
measures - see e.g. |CMS84j and [Sallbj and the references therein. The 
mere existence of such a do is evident for a compact manifold, however 
it can't be apriorily be supposed for a general manifold, hence has do be 
postulated. Moreover, to ensure a good approximation of curvature, this 
density factor has to be properly chosen (see, e.g. |SAZ12] ). thus tighter 
estimates for the mesh of the triangulation can be obtained from (llSp along 
with better curvature approximation. No less importantly, an adequate 
choice of the vertices of the triangulation, also ensures, via the thickness 
property, the non-degeneracy of the manifold (and of its curvature measures) 
-see [Sallbj . 

Proof 1 The theorem follows readily from Theorem 3 of |St76cj . Indeed, in 
the two dimensional case, the so called maximum and minimum curvatures, 
A;+, respective /c_ (see [St 76c] . p. 12, for the precise definitions) at the 
vertices of Mp^ coincide with the combinatorial Gauss curvature. Moreover, 
conditions (1) and (2) of Theorem 3 of jSt76c| are, due to the fact that 
here we are concerned solely with 2-dimensional simplicial complexes {PL 
manifolds), equivalent to our conditions (2) and (1) above, respectively. 
Therefore, the conditions in the statement of Theorem 3, |St76c] are satisfied 
and, by (ii) of the said result, the theorem above follows immediately. 

D 

Remark 3.3. It is easy to see that the theorem above extends to more general 
polyhedral surfaces. Indeed, by their very definition such surfaces admit 
simplicial subdivisions. However, during this subdivision, kj^, respective 
k- do not change, since the only relevant contributions to these quantities 
occur at the vertices, and depend only on the angles at these vertices, more 



l°in Stone's formulation ( [StTS] , p. 1062). 
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precisely on the normal geometry (see |St76c] . p. 12), that suffer no change 
during the subdivision process. 

Remark 3.4. The bound (jlSp is rather weak, as compared to the one for 
the classical case, but it is the only one supplied by Stone's result we made 
appeal to, namely Theorem 3 of jSt76c| . 

The proof above suffers from the disadvantage of making use of Stone's 
maximum and minimum curvatures (even though, in this context making 
appeal to them is rather natural). We can, however, provide a different 
proof, independent of Stone's work, but at the price of using some heavy 
(albeit classical) machinery, that, moreover, takes us away, so to say, from 
the discrete methods. (On the other hand, smooth, analytical tools are far 
more familiar to a large research community in CAGD, Imaging, etc.) 

Proof 2 The basic idea (which we first employed in |Salla) ) is to consider 
a sm^oothing M^ of Mpj^. Since, by |Mu66j . Theorem 4.8, smoothings ap- 
proximate arbitrarily well both distances and angled ^n on Mp^, defects are 
also arbitrarily well approximated. Given that the combinatorial curvature 
of Mp^ is bounded from below, it follows that so will be the sectional (i.e. 
Gauss) curvature of M^ . 

Unfortunately, the Gaussian curvature of M^ is positive only on isolated 
points (the set of vertices of Mj,^), so we can not apply the classical Bonnet 
theorem yet. However, we can ensure that M^ is arbitrarily close to a 
smooth surface -/Vf|., having curvature Gaussian curvature K{M-^) > Ol i 
Therefore, the classical Bonnet Theorem can be applied for M^, hence Mpj^ 
is compact and its diameter has the same upper bound (again using the 
same arguments as before 1) as that of M'^ (and M^), namely 

(17) diam(M|,i) < 



Kn 



More precisely, they are ^-approximation and, for 5 small enough, also e- 
approximations of Mpj^ - for details see [Mu66] . or, just for the minimal required facts, 
the Appendix of [Sal la] . 

This is easily seen by adding spherical "roofs" (of low curvature) over the faces, and 
then slightly modifying the construction, to ensure that the curvature will be positive also 
on the "sutures" of the said roofs, corresponding to the edges of the original PL manifold. 

i.e. S- and e-approximations 
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D 

Remark 3.5. Apparently, the bound for diameter given by the proof above, 
is tighter than the one obtained by Stone in [St76c| . However, we should 
keep in mind that, in practice, one is more likely to encounter PL surfaces 
as approximations of smooth oneso However, the larger the mesh of the 
approximating surface (i.e. the "rougher" the approximation), the larger 
the deviation of the approximating triangles from the tangent planes (at the 
vertices), hence the more likely is to obtain large combinatorial curvature. 
Hence, there is a correlation between size of the simplices and curvature, 
more precisely, the lower bounds in (jlSp . the lower ones in ()17p . 

Since the leitmotif of the previous section was metric (Wald) curvature, 
it is natural to ask whether a fitting version of the Bonnet-Meyers Theorem 
exists for this type of curvature? The answer is - at least in dimension 2 - 
positive: we can, indeed state an analogue of Meyers' Theorem, in terms of 
the Wald curvature: 

Theorem 3.6 (Bonnet-Meyers for PL 2-manifolds - Metric). Let Mp^ be 
a complete, 2-dim,ensional PL manifold without boundary, such that 

(V) There exists d^ > 0, such that mesh(Mp^) < di 

(ii') Kw{Ml^) >Ko>0. 

Then Mp^ is compact and, moreover 

(18) diam(M|.i) < 



0; 



Proof. We employ again the basic argument first used in |Salla| : Since 
distances (and angles) are arbitrarily well approximated by smoothings, it 
follows that so are metric quadruples (including their angles), hence so is 
Wald curvature. By [Bl53] (see also [BM70], Theorems 11.2 and 11.3), the 
Wald curvature at any point of non-trivial geometry M^, namely at a vertex 
V, K\y{v) equals the classical (Gauss) curvature K{v) (and, of course, this 
is also true a fortiori at all the other points, where both the smooth and the 
PL manifold are flat). Therefore the Gauss curvature of M^ approximates 
arbitrarily well the Wald curvature of Mp^, hence we can apply the same 
argument as in Proof 2 above to show that Mpj^ is, indeed, compact and, 
furthermore, satisfies the upper bound (fTTl) . D 



and, obviously, PL surfaces are PL approximations of their own smoothings 
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Remark 3.7. Like the previous theorem, the result above can be extended 
to polyhedral manifolds, and even in a more direct fashion, since Wald 
curvature does not take into account the number of sides of the faces incident 
to a vertex, but only their lengths. 

Remark 3.8. This result, as well as its generalization to higher dimensions 
(see l3.10() is hardly surprising, given the fact that, by |BGP92| . Theorem 3.6, 
Myers' theorem holds for general Alexandrov spaces of curvature > Kq > 0, 
and since Wald-Berestovskii curvature is essentially equivalent to the Rinow 
curvature (see [B153J ). hence to the Alexandrov curvature (see, e.g. |Gr99| . 
Chapter 1). Rather, we give, in the special case of PL surfaces (manifolds) 
a simpler, more intuitive proof of the Burago-Gromov-Perelman extension 
of Meyers' Theorem. 

In higher dimension, none of the arguments applied in both proofs of 
Theorem l3.1l are applicable, at least not without imposing further conditions: 

• Regarding the first proof: 

— In dimensions higher than 2, k^ and A;_ do not, necessarily equal 
each other (see |St76cj . Example 4, p. 14) and, a fortiori, they 
fail to equal the combinatorial Gauss curvature. They do, how- 
ever, according to Stone |St76c] . resemble in their behaviour the 
minimum, respective maximum sectional curvature at a point 
common to two 2-planes, that contain a given (fixed) tangent 
vector at the point in question. 

An important proviso should be added, however: While for the 
general PL simplicial complexes, the equality between A;+ and 
k- fails to hold, it is true for the most relevant - at least as 
far as our analysis is concerned - case of PL manifolds without 
boundary (see |St76cj . Example 3, p. 13). Consequently, it is 
not clear how to connect our proposed metric discretization of 
Ricci curvature with the the maximal and minimal curvatures of 
Stone (hence to combinatorial curvature, whenever they equal 
it - and each other)! ^ 



A natural attempt would be to use straightforward extensions of k+ and k- - let's 
denote them, for convenience, Ricmin and Ricmax. However, it is not clear (at least at this 
point in time) how expressive these definitions would prove to be. 
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Remark 3.9. It is true that the lower bound on kj^, as considered 
in Theorem 3 of |St76c] has a simple expression, in any dimen- 
sion, via a topological condition (cf. Lemma 5.1 of [St76cj ). 
namely that the intersection of any {PL) geodesic segment of 
ends p and q with the 2-skeleton of M|,^ is precisely the set 
{p, q} (with the exception, of course, of the case when the seg- 
ment is contained in a simplex (of Mp^). However, since the 
metric information contained in this new condition is void (or 
rather thoroughly encrypted, so to say) it has no apparent ad- 
vantage for application in conjunction with metric curvature. 

— For an application of the Stone's methods in conjunction with 
the metric curvature approach to any dimension, one would 
have to make appeal to Jacobi fields, as defined in |St76a] . How- 
ever, as discussed in the previous section, this would probably 
led to numerical instability. 

• As far as the second proof is concerned: 

— No smoothing of a PL manifold necessarily exists in dimension 
higher than 4 and, even if it exists, it is not necessarily unique 
(starting from dimension 4) - see |Mu60] . 

However, if such a smoothing exists, then the second proof of 
Theorem 13.11 (and of Theorem 13. 6p extends to any dimension, 
and we obtain the following PL (metric) versions of the classical 
results: 

Theorem 3.10 {PL Bonnet - metric). Let Mp^ he a complete, 
n-dmensional PL, sm,oothable m,anifold without boundary, such 
that 

(i') There exists do > 0, such that mesh(Mp^) < do; 
(ii')Kw{M^^)>K^>Q, 

where Kw{Mp^) denotes the sectional curvature of the "com- 
binatorial sections" i.e. the cells Ci (see Section 1 above). 
Then Mp^ is compact and, moreover 



(19) diam(M|,^) < 



vr 



In all honesty, we should add that the "rounding" argument 
of Proof 2 of Theorem 13.11 is not easy to extend directly — 
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if at all - to higher dimension. Instead, a more subtle argu- 
ment has to be devised. To this end we make appeal again 
to Stone's paper [St76c| . and we build the spherical simplicial 
complex M^ 1^ associated to the given PL (or rather piecewise- 
flat) complex Mp^. This is built as follows: Consider the 
sphere of radius R = R{cr) and radius O = 0{a), circum- 
scribed to a given simplex a, and its image a* = a*{R*) on 
a sphere of radius R* = R*{(7), R* > R, via the central projec- 
tion from O. Then we denote by simplicial complex obtained 
by remetrization of Mpj^ by the replacement of each a by its 
spherical counterpart a*. Then, by Lemma 5.5 of |St76c] . for 
large enough R* > R, the following holds for any pair of points 
p,q £ Mp^: distA/" ip,q) < CdistAf" ip*,Q*), for a certain 
constant C, where p*,q* denote the spherical images of p,q. 
Since the curvature at each vertex of the spherical simplex 
obtained by central projection of the simplices of Mpj^ onto 
their circumscribed spheres is smaller than the corresponding 
one (at the same vertex) in the PL (piecewise flat manifold), 
this holds a fortiori for M^^ . It follows from the classical 
Bonnet theorem (after applying the necessary smoothing) that 
diam(M]^^) < diam(M^p^^^). 

On the other hand, if we approach the problem of PL Ricci 
curvature from the viewpoint of the first part of the paper, that 
is of PL approximations of smooth manifolds, then the situa- 
tion changes dramatically. Indeed, even when such a smoothing 
M" (n > 3) exists, it is not probable that its sections provided 
in this manner by Mpj^ suffice to approximate well enough - 
let alone reconstruct - the Ricci curvature of M". In simple 
words, "there are not enough directions" in Mpj^ to allow us 
to infer from the metric curvatures of a PL approximation, 
those of a given smooth manifold M" (in fact, not not even a 
good approximation), hence we are faced again with a problem 
that we already mentioned in conjunction with the first proof, 
namely that of insufficient "sampling of directions" in PL ap- 
proximations. (On the other hand, increasing of the number of 
directions, i.e. of 2-dimensional sections (simplices) generates 
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a decrease of the the precision of the approximation, due to 
the (possible) loss of thickness of the triangulation - a problem 
which we have discussed in some detail in [Sallbj .) 

Remark 3.11. The considerations above show us that, unfortu- 
nately, no analogue in higher dimensions of the Myers' Theorem 
can be obtained by applying smoothing arguments). It is true 
that a Ricci curvature of the smooth manifold M" is obtained 
in terms of that of Mpj^, however, it is not clear, in view of the 
paucity of sectional directions (i.e. possible 2-sections), how 
precisely is this connected to its discrete counterpart. There- 
fore, we can obtain, at best, an approximation result (with lim- 
its imposed by the thickness constraint - see discussion above). 

We conclude with the following remarks: From the discussion above is 
transparent that, unfortunately, at this point in time, we can offer no proof 
for the general case, that is for non-smoothable PL manifolds of dimension 
n > A. To obtain such a proof for Bonnet's Theorem, one should adapt 
Stone's methods, as developed in [St76cj . while for a comprehensive gener- 
alization of Myers' theorem, one has the apparently more difficult task of 
adapting the purely combinatorial methods of |St76a] to the metric case. A 
quite different approach, but one that would allow us to extend the metric 
approach to quite general weighted CW complexes, would be to adapt For- 
man's methods developed in |Fo03j to our case. The essential step in this 
direction would be to find relevant geometric content (e.g. lengths, area, 
volume) for Forman's "standard weights" associated to each cell. 
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